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ABSTRACT 

Context. The Virgo direction has been observed at many wavelengths in the recent years, in particular in the ultraviolet with GALEX. 
The far ultraviolet (FUV) diffuse light detected by GALEX bears interesting information on the large scale distribution of Galactic 
dust, owing to the GALEX FUV band sensitivity and resolution. 

Aims. We aim to characterise the ultraviolet large scale distribution of diffuse emission in the Virgo direction. A map of this emission 
may become useful for various studies by identifying regions where dust affects observations by either scattering light or absorbing 
radiation. 

Methods. We construct mosaics of the FUV and near ultraviolet diffuse emission over a large sky region (RA 12 to 13 hours, DEC 0 
to 20 degrees) surrounding the Virgo cluster, using all the GALEX available data in the area. We test for the first time the utilisation 
of the FUV diffuse light as a Galactic extinction E(B - V) tracer. 

Results. The FUV diffuse light scattered on cirrus reveals details in their geometry. Despite a large dispersion, the FUV diffuse light 
correlates roughly with other Galactic dust tracers (coming from IRAS, Herschel, Planck ), offering an opportunity to use the FUV 
emission to locate them in future studies with a better resolution (about 5 arcsec native resolution, 20 arcsec pixels maps presented in 
this paper) than several usual tracers. Estimating the Galactic dust extinction on the basis of this emission allows us to find a smaller 
dispersion in the NUV-i colour of background galaxies at a given E(B-V) than with other tracers. The diffuse light mosaics obtained 
in this work are made publicly available.* 

Key words, dust, extinction, Ultraviolet: ISM,Ultraviolet: galaxies 


1. Introduction 


Large scale diffuse structures have been found early-on in the 
Galaxy (de Vaucouleurs| |1955| ), |Sandage| ( |1976[ ) publishing for 
the first time detailed optical images of high latitude clouds 
found in a Galactic polar caps survey. However, It is only with 
the IRAS satellite that we were able to discover a rich geometry 
of far infrared emitting cirrus clouds over the full sky (Low et 
al.|1984| ). Deep optical imaging then revealed their emission due 


to dust-scattering, showing a good correspondence with the far 
infrared (e.g.|de Vries & Le Poole|1985} |Guh athakurt a & Cutri 


|1994 ; Witt et al. 2008). The cirrus are also seen in the far ultra- 


* The FUV and NUV diffuse light mosaics will be available on the 
GUViCS web site http://galex.lam.fr/guvics/mosaics.html and at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 


violet (FUV) domain, in which a major component of the light 
results from dust scattering ( |Witt et al.|1997| ). 


Early studies of the FUV diffuse light 

can be found in Lil- 

lie & Witt 

(1976); Paresce et al. 

(1980);]Jakobsen et al. ( 

1984); 

Murthy et al. (1989 );'Fix et al. 

(1989); jHurwitz et al. ( 

1991); 

Perault et al. 

(1091): Murthy et al. (1999). 

Haikala et al. (1995 ) 


discovered with the FAUST telescope Galactic cirrus directly 
from their FUV emission. With the GALEX telescope, and its 
All sky Imaging Survey (AIS), exte nsive studies of the diffuse 
FUV emission became possible (e.g. Sujatha et al.|2010[ [Ham¬ 
den et al.|2013[|Murthy|2014| ). 


Studying the cirrus in FUV is useful for many reasons: 


UV data allows the determination of the scattering dust prop¬ 
erties in this wavelength domain, constraining dust grain 
properties (e.g. |Witt etal.l 19971 ISujatha et al.|2005) . 
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Rece ntly, Seo n et aL] ( |2014| ) and |Hod ges-Kluck & Bregman 
( |2014| ) have found ultraviolet haloes in edge-on galaxies that 
they interpret as reflection nebula from dust outside galax¬ 
ies. While it is a different component than the cirrus them¬ 
selves, characterising the dust properties in our backyard is 
also useful for such extra-galactic studies, provided the dust 
is of similar nature. 

Several works found an excess red emission (between 6000 
and 8000 A) in reflection nebulae, resulting from a lu¬ 
minescence excited by the FUV radiation of illuminating 
stars ([Witt & Schild] 1985|. [Guhathakurta & Tyson] (|1989|); 


Guhathakurta & Cutri (1994 ); Gordon et al. (1998); Szomoru 
& Gu hathakurta ( 1998)T |Witt et al. ( |2008 ) have shown that 
such Extended Red Emission (ERE) is also found in the dif¬ 
fuse interstellar medium and cirrus clouds. 

The presence of cirrus in our Galaxy directly affects the 
colour or detection of objects in deep observations of extra- 
galactic sources, since cirrus modify the background to 
which the objects are contrasted. It is especially true for low 
surface brightness regions in nearby galaxies such as tidal 
features that may be hard to distinguish from Galactic dust 
emission, or for lensing studies. This point is well illustrated 
by the Fig. 10 of [Due et aL| ( |2015| ). 

The presence of Galactic cirrus also affects UV/optical stud- 
ies of extra-gala ctic sources via dust extinction (e.g. |Bur stein 
& Heiles |1982|. Usually, data are simply corrected using 


Schle g eTet al.| ( T998| maps based on the 100 microns dust 
emission. 


Recently, a large area around the Virgo cluster has been sur¬ 
veyed at many wavelengths: especially by GALEX in the FUV 
and near ultraviolet (NUV) bands forming the GUViCS sur- 
vey ( [Boselli et al.||2011| ). These data have been used already 
for several works studying the effect of the cluster environment 
on the evolution of galaxies (Cortese et al. 2011| |Boselli et al 


2014 ), the effect of ram-pressure on star formation (Boissier et 


al.|2012| ), or stripped systems (Arrigoni Batta ia et al.|2012| ). The 


source catalogs have been publised in Voyer et al. (2014). 

The area has also been observed in the frameworks of the 
Next Generation Virgo Cluster Survey (NGVS) in the optical 
(Ferrares e et al.| 2012), of the Herschel Virgo Cluster Survey 
(HeViCS) in the far infrared ( [Davies et al.| 2010| ), and o f the 
Arecibo Legacy Fast ALFA survey (ALFALFA) in HI ( |Gio-| 
vanelli et al. 2005). Deep B and V images (reaching 28 .5 mag 


Mihos et al. 


(2015) 


arcsec - ^ in the B band) were obtained by 
in about 15 deg 2 around M87 and M89, using the Schmidt tele¬ 
scope at Kitt Peak. While several of these surveys are mostly 
dedicated to extra-galactic science, they can also be used to study 
the cirrus component in our Galaxy, or may be affected by them. 

In this paper, we concentrate on the UV data from the GU¬ 
ViCS survey. We characterise the cirrus in the UV and com- 
pareour findings to other wavelengths. Recently, Hamd en~et| 
[ah ( j2013j ) performed a full sky study of the ultraviolet diffuse 
light in the GALEX All sky Imaging Survey. We focus instead 
on the Virgo cluster area, taking advantage of the large num¬ 
ber of observations performed allowing us to use deeper data, 
and to visually inspect FUV detected cirrus and compare them 
to ancillary data in a small region of the sky, while their work 
is statistical in nature. This will be useful to any analysis based 
on deep data (from the many recent surveys of Virgo) that may 
be subject to pollution by diffuse low surface brightness emis¬ 


sion. For instance, our FUV maps were already used by Dur 
rell et al.| |2014| to make sure their results concerning the spa¬ 


tial distribution of globular clusters in Virgo were not affected 
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Fig. 1. Objects-subtracted NUV counts (based on 1 aremin pixels, 
median-filtered in 3x3 pixels boxes) as a function of the ecliptic lati¬ 
tude. The numbers on the right axis indicate the surface brightness in 
units of AB magnitude per square arcsec. The colorbar indicates the 
decimal logarithm of the number of pixels per 2D bin. The red line 
shows the average. The grey line is the adopted background (Sect. |3.2| ). 
The magenta squares report predictions of the GALEX zodiacal light 
tool. 


by cirrus. The presence of cirrus may be crucial in the study of 
low surface brightness features that may be associated to tidal 
interactions, expected to be frequent in the cluster environment 
of Virgo. A school case is the interacting pair NGC 4435/4438. 
[Cortese et al.| ( |2010| ) studied this system with a combination of 
multi-wavelength data including far infrared and far ultraviolet 
observations. These data helped them to elucidate the nature of 
a “plume” detected in deep optical data. The cirrus are also eas¬ 
ily recognized in a simple visual examination of the deep optical 
imaging of [Mihos et al.| ( [2015| ). While they seem to correspond 
well to the location where we detect them in FUV or infrared, 
a detailed analysis (distinguishing the cirrus vs faint emission 
linked to the galaxies, and comparing to other wavelengths) will 
be performed in the future. 

In Sect. [2] we present our data and how we constructed the 
mosaics (made publicly available with this paper) over a large 
sky area (RA 12 to 13 hours, DEC 0 to 20 degrees) surrounding 
the Virgo cluster. An analysis (UV properties of cirrus, compar¬ 
ison with other wavelengths) is performed in Sect. [3j and our 
conclusions are summarised in Sect. [4] 


2. The GUViCS mosaics 

2.1. FUV and NUV GALEX mosaics 


The GUViCS survey was presented in [Boselli et al.| ( 2011| . It 
combines a large numbers of tiles observed with GALEX in the 
FUV and NUV bands (with effective wavelengths respectively 
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Fig. 2. Objects-subtracted FUV counts (based on 1 arcmin pixels, me¬ 
dian filtered in 3x3 pixels boxes) as a function of the ecliptic latitude. 
The numbers on the right axis indicate the surface brightness in units 
of AB magnitude per square arcsec.The colorbar indicates the decimal 
logarithm of the number of pixels per 2D bin. The red line shows the av¬ 
erage. The grey line is the adopted background (Sect. |3.2| ). The magenta 
squares report predictions of the GALEX zodiacal light tool. 


1528 and 2271 A, and a spatial resolution of about 5 arcsecs). 
The minimal exposure time observed corresponds to the AIS 
survey (exposure times of typically 100 seconds, reaching a sur¬ 
face brightness of about 26 mag arcsec -2 ) but are often typical 
of the Medium Imaging Survey or Nearby Galaxy Survey (1500 
seconds exposure times, 28 mag arcsec -2 ). A few tiles have even 
much deeper data because they were observed for other purposes 
and are included in GUViCS. An update with all the tiles we in¬ 
clude in our work and several source catalogues can be found in 
Voyer et al. ( |2014| ). The circular ~ 1.2 degree diameter field of 
GALEX offers a good basis to map large-scale structures. How¬ 
ever, Galactic cirrus extend beyond this size. To circumvent this 
difficulty, we took advantage of having access to all the files 
needed to generate mosaics of the full area of the GUViCS sur¬ 
vey in both the FUV and NUV band of the GALEX telescope. 

We started from the following GALEX pipeline products: in¬ 
tensity maps, exposure maps, object detection masks. After some 
tests, we decided to truncate them at a radius of 0.58 degrees to 
avoid problems often observed on the edge of GALEX fields 
(i.e. increased noise, elongated point-spread function). We also 
removed a few tiles originally in the survey but with obvious 
positioning issues (see Voyer et al. 2014[ for examples). 

As a first step, in all the tiles, we replaced pixels in the object 
detection maps by surrounding values. Then, we used the Mon¬ 
tage softwarqj to co-add the objects subtracted intensity maps 
by weighting each image by its exposure time. Montage allows 


http://montage.ipac.caltech.edu/ 


to project the final image on an arbitrary pixel grid. We tested 
several resolutions looking for a balance between resolution and 
signal to noise. We present in this paper maps a mosaic with 1 
arcmin size pixels over the full Virgo area for tiles from the AIS 
(having the advantage of being homogeneous and covering the 
largest possible area). We also produce an image with pixels of 
20 arcsecs size showing finer details in the background structure 
owing to deeper exposures (i.e. exposure times of at least the 
1500 seconds typical of the GALEX Medium Imaging Survey). 

This method is similar to the one used in lHamden et al.l 
P013| ). The main differences are : they rely on sky background 
images from the GALEX pipeline while we replace the detected 
objects within the masks by neighbouring pixel values by our¬ 
selves (the pipeline background is smoothed and may not re¬ 
produce small-scale variations in the background). Second: they 
grow the masks of detected objects to limit the pollution of data 
by actual objects, while we keep the masks at their size to in¬ 
crease our chance to see real structures (even if the absolute flux 
level may be in some case polluted). Finally, their final map has 
pixels covering 11.79 arcmin 2 , while we choose 1 arcmin 2 and 
20x20 arcsec 2 pixels in our study of a more limited area. The 
higher resolution allows us to see details of the structures re¬ 
vealed in the FUV with deeper data than the All Sky Imaging 
Survey that they use in their full sky statistical study. 

The FUV mosaics obtained in this way presented a few arte¬ 
facts due to the absence of data in a small area on the edge of the 
FUV detector. We edited manually the FUV mosaics to inter¬ 
polate over such regions. During this visual inspection, we also 
manually masked a few sources that had been badly subtracted 
(usually very bright stars or extended galaxies). These correc¬ 
tions concern a small number of pixels. We checked that they 
do not affect the conclusions presented here (we did perform our 
analysis with both the raw mosaic and the edited image). The 
only visible difference is a tail of FUV bright pixels correspond¬ 
ing to badly subtracted objects in the uncorrected map. 

As a final step, to obtain smoother maps, we found useful to 
median-filter our mosaics in 3 x 3 pixel boxes. This alleviates 
potential difficulties linked to the choice of the pixel grid origin 
and it reduces the noise at the very low surface brightness of 
the cirrus, at the price of degrading further the resolution of our 
mosaic to about 3 arcmin (for the lower resolution mosaic). This 
is the same order of magnitude as the resolution of other dust 
tracers to which we will compare our map. 

2.2. Units 

Different communities express the UV surface brightness in dif¬ 
ferent units (see |Leinert et al.|||1998| for an exhaustive list of 
surface brightness units and their conversions). To make it easy 
to convert between various units used in this paper, we express 
most of our maps in GALEX counts (per sec). To have gras- 
pable numbers, we used counts per arcmin 2 in most cases. When 
FUV/NUV comparison was involved, we switched to AB mag¬ 
nitudes and surface brightnesses which are relatively standard. 
Conversion between counts and surface brightness is straight¬ 
forward 

FUV AB magnitude = -2.5/<9gl0(counts sec -1 ) + 18.82 (1) 


NUV AB magnitude = -2.5/<9gl0(counts sec l ) + 20.08 (2) 

2 http://galexgi.gsfc.nasa.gov/docs/galex/FAQ/counts_ 
background.html 
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Fig. 3. Left: mosaic showing the diffuse FUV emission obtained from AIS data (1 arcmin pixels, median smoothed on 3x3 pixels). Right: diffuse 
FUY emission obtained on the basis of deeper da ta (o n 20 arcsec pixels, median smoothed on 3x3 pixels). The contours indicate roughly the 
regions in which the cirrus emission is found (Sect. |3.2| ) 


"Continuum Units" (CUs, photons cm -2 s -1 sr -1 A -1 ) are of¬ 
ten used in studies of the FUV diffuse light. It is straightforward 
to convert from counts to CUs for the FUV GALEX band: 

counts sec - Arcmin -1 = CUs/1260 (3) 


2.3. Other wavelength maps 


Large scale variations in the sky background are known to be 
related to zodiacal light, and cirrus. These components are often 
studied by looking at their emission in the far infrared. |Schlegel| 
|et al.| (19^8] ) produced maps of IRAS/COBE 100 jum emission 
that have been widely used in the literature. The IRAS data were 
reprocessed by |Miville-Deschenes & Lagache| ( ]2005| ) to obtain 
at 100 microns a slightly improved angular resolution (4.3 ar¬ 
cmin), and a better calibration, taking into account the time and 
brightness response of the detector. This improved map is called 
IRIS (Improved Reprocessing of the IRAS Survey). [Planck CoL] 
|laboration et al.| ( |20l3] ), however, mention that the zodiacal light 
subtraction was done differently in the IRIS project with respect 
to |Schlegel et al.| ( |1998| ), and proposed a simple correction to 
use the large scale variations from Schlegel et al.|(|1998|) and 


the small-scale variations from 


Miville-Deschenes & Lagache 


( 2005| ). We use the IRIS 100 microns map and adopt the same 


correction in our study. Finally, we prepared a mask with all 
VCC galaxies detected by IRAS and with some other bright ob¬ 
jects such as known field galaxies and stars. The mask excludes 
up to twice the 100 microns IRAS resolution around each source. 
Masked regions were interpolated over. 

The Planck collaboration has used IRAS and Planck all sky 
maps to compute the properties of the Galactic dust thermal 
emission. The collaboration made available the dust maps, in¬ 


cluding dust opacity, temperature and radiance (Planck Collab- 
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|oration et al.|2011[|2013| ). The dust radiance was converted into 
a reddening E{B - V ) ma£] using the colours of SDSS quasars 
as a reference. This map should be adequate to estimate the dust 
attenuation for extra-galactic sources, with a spatial resolution of 
about 5 arcmin. 


The HeViCS project (Davies et al.j 2010| collected Herschel 
data in 4 fields within the area studied in this paper. In this work, 
we use the 250yum maps observed with the SPIRE instrument 
(spatial resolution of about 18 arcsecs). The data was reduced us¬ 
ing the Herschel dedicated software HIPE, adopting the vl 1 cal¬ 
ibration tree, assuming the gain correction for extended sources, 
and combining the scans with the built-in destriper. This was 
found to preserve the large scale diffuse emission better than 
the dedicated pipeline used in other HeViCS works, which was 
tailored for extra-galactic targets (see, e.g. Auld et al. 2013). 
Further details on the data reduction, together with a study of 
the properties of the diffuse cirrus emission will be presented in 
Bianchi et al. (in preparation). We interpolated these maps over 
the same mask as the one used for the IRAS 100 microns emis¬ 
sion. However after this step, the HeViCS tiles are still filled 
with many faint sources. The background was thus determined 
following the same procedure as in Pappalard o et al. | ( [2014| ), i.e. 


with SExtractor. The cirrus 
Baes et al. ( |2014| ) discussed 


computing iteratively a backgrounc 
are well visible in HeViCS images 
their role as potential non extra Galactic sources in their cross¬ 
correlation of HeViCS and Planck maps. A HeViCS study of the 
cirrus is in preparation (Bianchi et al.), thus the HeViCS find¬ 
ings are not discussed here in details but the data are shown for 
illustrative and comparison purpose. 


3 http://irsa.ipac.caltech.edu/data/Planck/release_l/ 
all-sky-maps/ 
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Fig. 4. Gallery of multi-wavelength tracers of the cirrus in the Virgo Area. Top-Left: composition of FUV (blue) and NUV (red) mosaics from 
GUViCS. Bottom-Left: Planck E(B - V ) map . Top-Right: 100 microns IRIS map. Bottom-Right: 250 microns HeViCS map. The yellow contours 
are the broadly defined regions of cirrus detected in FUV (Sect. |3.2| ). 


We extracted the parts of the 100 microns and E(B - V ) 
Planck map covering our UV mosaics using the Aladin sky atlas 
(Bonnarel et al. 2000| ). The 250 micron maps directly come from 
the HeViCS collaboration. All were registered to the same pixel 
scale as our various UV maps (1 arcmin and 20 arcsec). 


In principle, it would be very interesting to compare the ul¬ 
traviolet detection of cirrus to the all sky dust emission computed 
in the WISE W3 band by |Meisner & Finkbeiner| ( [2014| ). At 12 
microns, this map reveals the emission from small grains and 
Polycyclic Aromatic Hydrocarbon (PAHs) absorbing mostly ra¬ 
diation in the ultraviolet. Scattered ultraviolet light should then 
spatially occur close to the 12 microns emission (while the ther¬ 


mal dust is heated by the general interstellar radiation field that 
can be spatially mismatched from the ultraviolet strong sources). 
We did download the dust emission maps published by Meisner 
|& Fi nkbeiner ( |2014| ). Sadly, the Virgo area is crossed by a wide 
stripe of corrupted data, making them unusable for a statistical 
comparison that would correspond to our goals. A detailed com¬ 
parison of some regions could still be possible in the future (we 
illustrate this point in section]?]). 
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3. Analysis of large scale variations 

3.1. Zodiacal light 


For general studies of the zodiacal light in the UV, the reader is 
invited to consult e.g. |Frey et al.| ( pW7] ); |Murthy| ( |2Q14| ). Here 
we only analyse its effects on our mosaics and their contamina¬ 
tion on the cirrus emission. In the NUV band, we found that our 
object-subtracted map presented a strong gradient in its back¬ 
ground, which can be ascribed to the zodiacal light. The GALEX 
site propo ses a zodiacal light tooQ based on the zodiacal light 
study by |Leinert et al.| ( |1998| ). It predicts the level of zodiacal 
light expected for a position and an epoch. We checked with this 
tool that indeed a global decreasing gradient is expected in the 
NUV observations along the ecliptic latitude. 

We show in Fig. [T] the actual gradient with ecliptic latitude 
that we measured in our GALEX diffuse NUV image over the 
full Virgo cluster area. It is compared to the prediction of the 
GALEX zodiacal tool for a few coordinates regularly covering 
the area, adopting the date of January 1st. Variations with the 
date do occur, but change in the same way at all latitudes. These 
variations are at most time lower than about 25 percent, except 
during a zodiacal emission spike, when GALEX observations 
were not executed. At low ecliptic latitude, we find large val¬ 
ues of the background, consistent with the GALEX tool model. 
At larger latitude, our background decreases, as the model does, 
but by a smaller amount. This may be due to either an under¬ 
estimation of the zodiacal light model, or that other sources start 
to contribute to the background (ghosts of stars, artefacts,...). We 
note, however, that if the background was dominated by other 
sources, it should not depend on the latitude. The grey line in 
Fig. [T] shows the gradient that we determined by hand to approx¬ 
imative^ pass by the darkest regions within the map at different 
latitude. This gradient is only used in the next section to consider 
the cirrus NUV detection and in Fig. [4] 

We proceeded to the same exercise in FUV (Fig. In this 
case, the expected level of the zodiacal light is extremely low. 
The observed level is much higher than the prediction and is 
independent of the latitude (as found by Hamden et al. 2013 


on a larger scale). Other effects thus dominate the backgrounc 
(e.g. ghosts, residuals from subtracted objects, unresolved dis¬ 
tant sources, and Galactic cirrus studied in the following). 

We did also look at how the diffuse FUV and NUV back¬ 
ground counts vary with the Galactic latitude. We found no trend 
with the Galactic latitude in FUV, but we investigate a limited 
range of Galactic latitude, above 55 deg, while most of the trend 
due to cirrus is expected to occur at lower latitudes ( Ham den et| 
|ak]|2013| ). An anti-correlation is found in the NUV background 
counts, but we interpret it as the result of the variation with the 
ecliptic latitude (the zodiacal light described above), since in our 
field, the Galactic and ecliptic latitudes vary in the same direc¬ 
tion. 


3.2. Ultraviolet cirrus detection 

In the FUV, we measured the background level in several re¬ 
gions among the darkest area of the sky in our mosaic. We did 
this by using the IRAF task IMEXAMINE in several rectangu¬ 
lar regions, in a similar manner as previous works with GALEX 


data (e.g. 

Boissier et al. 2007, 2008 

2012). As it can be guessed 

from Fig. 

2j this background is found 

at a level of 0.25 counts per 


sec per arcmin 2 (with a 0.11 3-cr uncertainty per pixel). We cre¬ 
ated a FUV backgroud-subtracted map by subtracting this con- 

4 http://sherpa.caltech.edu/gips/tools/zlct.html 
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stant value to the FUV mosaic. By proceeding this way, we ba¬ 
sically remove any emission not related to the cirrus themselves, 
that could be due to air-glow or extra-galactic sources, and that 
should correspond to the "offset" observed by [Hamden et al. 
( |2013[ ). We acknowledge that this background could include a 
very diffuse cirrus contribution that is then subtracted in our ap¬ 
proach. The nature of this diffuse component emission is still 
unknown and the subject of several studies (see e.g. |Henry et al. 
|2015[ who test the hypothesis of a dark-matter interaction with 
the interstellar medium, or of small dust grains; and references 
therein). 

The remaining FUV emission above this constant back¬ 
ground is assumed to be due to cirrus. It is true that some part of 
it could have other sources (extended halos of galaxies, tails or 
streams, intracluster light) but the result discussed in the rest of 
the paper seems to be consistent with this emission being dom¬ 
inated by the Galactic cirrus. We created a detection map by 
adopting a 3-cr threshold above this level. This map was then 
smoothed (by a Gaussian filter of 6 arcsec) with the only goal 
to compute contours (in ds9) at a single level (0.1 counts sec -1 
arcmin -2 ) in order to geometrically identify broadly the cirrus 
regions. These contours are not directly used in the science anal¬ 
ysis of this paper, but we overlay them in our figures. They are 
shown in Fig. [3] over the the FUV mosaics obtained from AIS 
data and from deeper GUViCS data. 

In the NUV, we first constructed a background map to re¬ 
produce the gradient observed in Fig. |T] We used the darkest re¬ 
gions to estimate its zero-point. This background was subtracted 
to the NUV image. We then inspected the sky level within this 
background-subtracted image, and found it to be flat. Measuring 
the sky level in dark regions we obtained the value of 0.1 counts 
per sec per arcmin 2 with a 3-cr value of 0.52: it is consistent with 
zero, demonstrating that the sky gradient subtraction procedure 
performed well. We then produced a NUV "detection map", by 
adopting also here a threshold at this 3-cr value. 


3.3. UV Properties of Galactic cirrus 

Fig [4] (top-left) shows the backgroud-subtracted images in NUV 
and FUV with the FUV detection contours. The comparison 
definitively shows the blue colour of the cirrus with respect to 
the background. To characterize the UV properties of our cirrus, 
we show in Fig.[5]the relation between the FUV-NUV colour and 
the FUV surface brightness, obtained from the object-subtracted, 
backgroud-subtracted 3x3 median-filtered 1 arcmin size pixel 
maps. The higher resolution maps (having a higher signal to 
noise, but smaller and inhomogeneous coverage) produce very 
similar results from this point of view. Our data point to bluer 
colours for brighter surface brightnesses. |Sujatha et al.| ( |2010| ) 
found such an increase in the FUV/NUV ratio that they inter¬ 
preted as a sign of fluorescent emission of molecular hydrogen. 
In our data, however, we cannot rule out a detection issue, as 
can be seen from Fig. [5] Taking into account the thresholds dis¬ 
cussed in Sect. |3.2| and restricting ourselves to the pixels where 
the signal to noise is larger than 3 at both wavelengths, we find 
that the typical FUV-NUV colour of the cirrus is close to 0. This 
is very similar to the colour of FUV-NUV ~ 0.11 that can be de¬ 
duced from the spectral distribution found by Gondhalekar et al. 
( 1980[ ) for the interstellar radiation field, probably typical of B- 
stars illuminating the dust clouds. In the scattered haloes of ex¬ 
ternal star-forming galaxies, |Hodges-Kluck & Bregman] ( |2014| ) 
also found quite blue colours with SEDs ( AF ^) that are flat (or 
declining from FUV to NUV), corresponding to a FUV-NUV 
colour of ~ 0.4 (or bluer). 
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Fig. 6. Comparison of several cirrus indicators in a zoomed region with deep FUV data (top-left). This map is compared with the 100 microns 
IRIS map (top-right), HeViCS 250 microns (bottom-right), Planck E(B - V) map (bottom-left). 


3.4. Comparison to other wavelengths 


Fig. [4] shows the 100 microns, 250 microns, and Planck dust red¬ 
dening maps compared to a FUV+NUV combination (obtained 
from the shallow data covering a large fraction of the full area). 
Globally, a similar morphology is found with the various trac¬ 
ers. Such an agreement between various indicators for the cirrus 
presence has been already discussed before e.g. by Guhatha kurta| 
& Cutri ( 1994|) for the optical, 100 microns and gas distribution, 


or 


Seon et al. (201 1) for relations between FUV, gas column den¬ 


sity, 100 microns radiation. The agreement between the various 
tracers is, however, not perfect. For instance, some regions of 
the Planck E(B - V ) or far infrared maps are dark whithin UV 
contours in Fig. [4] (e.g. around 12h55m, 9 degrees). Small-scale 
variations do bring clues about the physical properties of the dust 
( [Guhathakurta & Cutri||1994| even if another source of disper¬ 
sion is due to the anisotropic interstellar radiation field and the 
precise geo metry of dust and bright stars (Wi tt et al.|1997| |Seon| 
|et al.|20lT] ). It is beyond the scope of this paper to explore these 
aspects in details here, but the maps that we distribute may help 
to explore such studies in the future. 


We show in Fig. [6] a comparison of several dust tracers in a 
region for which we do have data deeper than the AIS, using our 
20 arcsec pixel map. It reveals that the UV allows to follow in 
details the structure of the cirrus with a better resolution than the 
traditional tracers. 


Beyond a simple visual comparison of the maps, we now 
proceed to a quantitative comparison, pixel by pixel, of the FUV 
cirrus emission and other dust tracers. Fig. [7] and [8] show the 
distributions of pixels, and Table 1 provides the correlation co¬ 
efficients between the FUV and other dust tracers pixels. 

A relation between FUV and 100 microns emission had been 


reported in several studies fPerault et al. 1991 Witt et al. 1997; 

Sasseen & Deharveng 

1996: 

I. We show the relation we find in 

the Virgo direction in Fig. [7] 

The relation found by Hamden et 


with a 3x3 pixels median-filter (an adequate pixel size taking into 
account the 100 microns resolution). In this comparison, we used 
only the slope of the 100 microns-FUV relation of |Hamden et al.| 
( [20131 ) since the intersect is believed to be due to uniform sources 
that we subtracted in our map construction. We adopt their slopes 
in the Galactic latitude ranges corresponding to Virgo (60-75 and 
75-90). The good agreement suggests that our background mea¬ 
surement and subtraction was efficient. On the other hand, we 
note an important dispersion, a sign that the dust emission in 
the far infrared and the FUV reflection on dust are non perfectly 
coincident. We found a similar relation between the FUV and 
the 250 microns emission, albeit with more dispersion. The 250 
map, however, covers only a part of our FUV data. 

The dispersion in the relation between the FUV and infrared 
pixel values (reflected by the relatively poor correlation coeffi¬ 
cients) is likely due to the fact that while the two tracers globally 
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Fig. 5. FUV-NUV colour as a function of surface brightness (based on 
the 1 arcmin pixel size mosaics). The lines shows our detection thresh¬ 
olds in FUV and NUV (shaded area are below the detection threshold 
and thus subject to large uncertainties). The colour-bar indicates the 
decimal logarithm of the number of pixels in each 2D bin. 


Table 1. Pearson correlation coefficients 


X 

Y 

Correlation coefficient 

variable 

variable 

(Pearson) 

FUV 

100 microns 

0.57 

FUV 

250 microns 

0.47 

FUV 

E(B - V ) Planck 

0.61 


trace the shape of the cirrus, they are more complementary than 
equivalent. This can be better seen in Fig.[9]in which we show the 
ratio of the FUV diffuse emission to the IRIS 100 microns emis¬ 
sion. We first applied a 6 arcmin Gaussian convolution to each 
image to reveal the differences on large scales and avoid being 
dominated by small scale fluctuations (or resolution effects). It 
is evident from this figure that this ratio is not identical in the 
field but varies with position. It is not only due to noise since it 
varies also within the contours indicating regions where FUV is 
definitively detected. Fig. [TO] shows a zoom in a region where a 
clear elongated cirrus structure is seen, part of it being obvious 
through its FUV emission, part of it through its far infrared emis¬ 
sion. Some of the mismatch could be due to the relative geometry 
of stars and various type of dust. The FUV scattered light may 
occur on small grains predominantly in regions close to FUV 
emitting stars, while the far infrared radiation corresponds to the 
thermal emission of the dust heated by the interstellar radiation 
field, heating that can also occur in the absence of the hottest 
stars in the immediate vicinity. 
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Fig. 7. FUV counts vs 100 microns intensity diagram (1100). The colour 
scale shows the logarithm of the numbers of pixels within 2D-bins. The 
grey lines indicate the slope found by |Hamden et al.| < |2013| ) for two 
Galactic latitudes corresponding to the GUViCS area. The maps used 
have 1 arcmin size pixels and have been median-filtered over 3x3 pixels. 
The 100 microns correspond to IRIS repro cessed IRAS data corrected 
on large scales as explained in section \23\ 


3.5. FUV as a dust tracer in the Virgo Cluster area 

We show the distribution of pixels in the FUV emission vs 
Planck Collaboration et al. ( [2013 1 ) E(B - V) extinction plane in 
Fig. [8] The two variables are correlated. A very similar trend is 
obtained on the basis of the Schlege fet al.| ( |1998] ) maps (the fig¬ 
ure includes the best linear fit in this case). This indicates that 
the FUV diffuse emission is at least loosely linked to the extinc¬ 
tion as measured through the dust emission. We explore quanti¬ 
tatively the link between the FUV diffuse emission and the ex¬ 
tinction in this subsection. The FUV map used is the 1 arcmin 
pixel one, median-smoothed over 3x3 pixels. The Planck extinc¬ 
tion map has a resolution of a few arcmin, thus this scale is perti¬ 
nent for our comparison (we however computed the distribution 
at higher resolution -with the 20 arcsec pixels map- and found 
very similar trends). As expected, since the FUV and the Planck 
maps do not correspond perfectly to each other, some scatter is 
found. The figure, however, suggests a clear average relation be¬ 
tween the two quantities. We computed averages within bins of 
FUV emission and performed a simple fit to these values to find: 


E(B -V) = 0.02378 + 0.1105 x FUV(counts s 1 arcmin z ) (4) 


- 2 \ 


A shallower slope (~ 0.09 instead of 0.11) is found if the 
Schlegel et al 7|([1998| ) map is used rather than the Planck Collab 


oration et al. 


( 2013| ) one. Considering the dispersion, the value 


of this slope is uncertain, but a correlation is clearly present. As¬ 
suming the underlying average relation between FUV emission 
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Fig. 8. Planck E(B - V ) vs FUV counts. The colour bar indicates the 
logarithm of the number of pixels within each 2D-bins. The dashed line 
shows averages within FUV flux bins. The solid line indicates a simple 
linear fit to the average trend. The dotted line indicates the result of a 
similar analysis based on the |Schlegel et al.| ( |1998j) E(B - V) map. 
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Fig. 9. Ratio of the FUV diffuse emission (in the 1 arcmin pixel map) 
to the 100 microns emission, after applying to both image a Gaussian 
filter of 6 arcmin. The contours indicate roughly the regions in which 
the cirrus emission is found (Sect. |3.2| ). 



Fig. 10. Combination the FUV emission (blue) and 100 microns emis¬ 
sion (red) along an extended cirrus feature, spreading from left to the 
right part of the figure. The scaling was chosen to show the comple¬ 
mentarity of the two dust tracers, since the colour is changing along the 
structure. 


and dust extinction holds at small scales, we invite researcher 
working in the Virgo area to use our distributed 20 arcsec (or 
1 arcmin) map and this relation to predict "alternative" extinc¬ 
tion values with respect to the traditional |Schlegel et al.| ( |1998| ) 
ones (or Planck ones). The advantage is that these FUV maps 
reveal detailed structures of the cirrus not visible in the long- 
wavelength maps (l ack of resolution/sensitivity, and comple¬ 
mentarity : see Sect. |3.4| ). We warn the reader however: a 1 to 
1 relation is not expected as the the FUV emission is sensitive 
on the geometry of dust and UV-emitting stars ( |Witt et al.|1997[ 
|Seon et al.|2011| ), Equation]?] should then be used with caution, 
especially for individual objects. Despite this drawback, we be¬ 
lieve that an attenuation computed from our FUV diffuse maps 
provides an alternative value of the reddening that may be com¬ 
plementary to traditional dust tracers. 


3.6. Testing the FUV extinction determination with an 
extra-galactic catalogue 

To test this new extinction estimation, we use background galax¬ 


ies detected in NUV using the GUViCS catalogue from Voyer et 
|aL| ( |2014| ). The idea to use background objects to study Galactic 
extinction is not new ( jBurstein & Heiles|1978| and was applied 


by Fukugita et al. |2004 t o SPS S data. In fact, the E(B - V ) 

| Planck Collaboration et al.|( |2013| ) map is based on the same idea 


applied to the colours of background quasars used to calibrate 
the relation between E(B - V) and the thermal dust emission. 
|Cuillandre et al.| ( [2001| ) applied a similar method to measure at¬ 
tenuation in M31 from background galaxies. They also explored 
the effect of attenuation on their colour and found it to be hard to 
detect. In our case, however, the used colour is more favourable 
as we expect a larger effect, using the NUV band, very sensitive 
to extinction. 

We do not determine the extinction from the galaxy counts, 
but rather verify that the colours of background galaxies are sta¬ 
tistically in agreement with the Galaxy dust reddening at their 
position. In Fig. [TTJ we show the trend between the NUV - i 
colour and reddening for Virgo and background galaxies, with 
a signal to noise larger than 10 in the Voyer et al. catalogue, 
leaving us with around 25000 galaxies (similar results are ob¬ 
tained when the minimal signal to noise is reduced to 2, with 
about 70000 objects). This number is of course much larger than 
the number of quasars in the same region. The galaxy colour 
is directly taken from the catalogue. The figure shows the av¬ 
erage and the dispersion of the colour of galaxies in various 
bins of Galactic reddening. We show four cases depending on 
which reddening is used: the one from the |Schlegel et al.| ( |l998| 
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maps, given in the catalogue; the one from the Planck E(B - V ) 
map; or the one derived from the FUV maps (the 1 arcmin 
and 20 arcsec pixel ones) applying equation [?| In all cases, the 
average colour becomes redder (with a broad dispersion) with 
the reddening. At least statistically, using the E(B - V) derived 
from the FUV produces the expected trend on the colour, based 
on a Milky Way extinction curve (|Cardelli et al.||1989|) with : 
A(NUV) = 8.74 E(B - V ) and A{i) = 1.98 E(B - V). It is hard 
to distinguish the relative performances of the various redden¬ 
ing maps since the dispersion in the colour of galaxies (linked 
to their star formation histories, metallicities, ages,redshift) is 
much larger than the reddening trend itself. The observed varia¬ 
tion of the NUV- i colour with E(B - V) is consistent with the 
one predicted simply with this Milky Way extinction law, sug¬ 
gesting that the various extinction tracers provide a similar level 
of accuracy. 

Looking into the details, however, it can be seen in Fig. |TT] 
that the dispersion of the colour at a given E(B - V) depends 
slightly on the source of dust reddening considered. The one de¬ 
rived from the FUV 20 arcsec map is smaller than the other ones 
by about 10 percent. The dispersion of the NUV - i colour of 
background galaxies in this diagram has two components, one is 
the intrinsic spread of colours of galaxies (due to their various 
star formation histories, metallicity, redshift) that is the same in 
each cases. The other one is obviously the incorrect assignment 
of the Galactic attenuation, which is changing depending on the 
dust tracer chosen. The “best” E(B-V) estimator should produce 
the smallest dispersion in this Figure (for a perfect estimator, we 
would be left only with the dispersion intrinsic to galaxies). We 
can thus conclude that the E(B - V ) derived from the FUV 20 
arcsec map does statistically a better job for our sample than 
the other ones, although the difference is small. We verified that 
this result is not due to the smaller area covered by this mo¬ 
saic: repeating the exercise on the 1 arcmin map restricted to the 
area used in the 20 arcsec study does not reduce its dispersion. 
We thus have to attribute this better result to the resolution and 
depth gain when using the 20 arcsec pixel FUV map to estimate 
E(B-V). 

Some caution should be used when using this attenuation es¬ 
timator. While we found that the diffuse FUV light correlates 
with E(B-V), and we use this fact to "calibrate" the redenning, 
it is clear that the FUV scattered light does not provide informa¬ 
tion on small grain absorption in FUV. Moreover, recent works 
suggest a large variability in the dust properties at high Galactic 
latitudes ( |Ysard et al.|2015| . [Peek & Schiminovich| ( |2013| ) hint 
in addition to deviations from the standard extinction curve in 
the FUV band. Deriving extinctions in the FUV band from our 
estimate may thus produce erroneous results. 


4. Conclusion 


We characterised the diffuse FUV emission in the Virgo clus¬ 
ter direction with the GUViCS data. We compared this emission 
(geometrically and statistically) to the NUV emission, IRIS 100 
microns emission, HeViCS 250 microns emission, E(B-V) from 
[Schlegel et al.|| ( [1998| ) and from the Planck collaboration. 

Cirrus are clearly detected in the FUV channel. They are 
characterised by a constant blue FUV-NUV ~ 0 colour (an ob¬ 
served trend with FUV surface brightness is in fact related to 
observational limits). 

We find an overall agreement between the locations of cir¬ 
rus seen in FUV with those from other tracers, although some 
distinctions are found; certain regions are visible in the FUV 
maps but do not appear obviously in the infra-red or Planck 


Article number, page 10 of[TT| 



Fig. 11. Top: NUV-i colour of SDSS galaxies in the Voyer et aL ([2014) 
catalogue as a function of the Galactic dust reddening (taking only ob¬ 
jects with signal to noise larger than 10). The solid lines show the av¬ 
erage, while the dashed lines indicate the 1-sigma dispersion. The dust 
reddening is derived either from the Schlegel maps (as given in the GU¬ 
ViCS catalogue), from the Planck E(B - V) map, or derived from our 
FUV 1 arcmin or 20 arcsec resolution diffuse emission maps. The ar¬ 
row indicates the expected reddening trend (see Sect. |3.6| ). Bottom: the 
1 sigma dispersion of the NUV-i colour is indicated as a function of 
the reddening. The range plotted in both panel corresponds to the bins 
where we do obtain E(B - V) values for each method. 


maps and vice versa. Despite these regions, a correlation is found 
when the various cirrus indicators are compared on a pixel basis 
(with, however, a large dispersion). The FUV reveals details in 
the structure of the cirrus that are not always seen at any other 
wavelength, owing to a good spatial resolution and sensitivity, 
but also to a complementary between the FUV and far infrared 
emission. 

The diffuse FUV emission maps produced in this work may 
be useful for many purpose: 

- The study of the dust properties in the UV (albedo, scat¬ 
tering), or the properties of cirrus in general with multi¬ 
wavelength maps matched at the same pixel size. For in¬ 
stance, a small region with both deep FUV data and good 
quality WISE W3 dust emission is shown in Fig. [12] The 
good agreement between the two tracers could confirm that 
the 12 microns emission comes from small grains and PAHs, 
which absorb the FUV radiation. 

- In conjunction with deep optical data, the study of the Ex¬ 
tended Red Emission phenomenon. 

- Identification of regions "polluted" by faint structures related 
to Galactic cirrus. This could be very useful for the exploita¬ 
tion of the surveys coming to an achievement in the Virgo 
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Fig. 12. A region around coordinates RA=12:36:07.201, DEC=+14:16:39.98 is shown as seen in the 20 arcsec p ixels FUV di ffuse light map (left), 
in the NGVS ( [Ferrarese et al.|2012| > g band image (middle), and in the WISE W3 diffuse emission (right) from Meisner & Finkbeiner (2014) on 
the same scale. To ease the comparison, a single contour was computed in FUV and overlaid on the other images. 


area, especially NGVS in the optical and GUViCS in the UV. 
We further illustrate this point in Fig. [12] showing a compari¬ 
son of the same area in our FUV diffuse map and in a deep g- 
band NGVS image (the brightest cirrus feature here reaches 
27.2 mag arcsec -2 ). This figure shows that the cirrus are also 
detected in these deep high quality optical images, owing to 
an efficient restoration of the true sky background of Mega- 
Cam data. The similar distribution of the structures at both 
wavelengths is convincing that it is well related to cirrus. 
Provide alternative extinction maps by using the correlation 
found between the diffuse FUV emission (as measured in the 


map we deliver) and the reddening, as detailed in section 3.5 


The 20 arcsec FUV diffuse emission is slightly better than 
other dust tracers based on the relation between E(B -V) and 
the NUV - i colour of background galaxies (section [U6| ). 


For this reason, and while there is certainly a lot of analy¬ 
sis and understanding to be still performed, we put these FUV 
maps (1 arcmin pixels based on the AIS, and 20 arcsec pix¬ 
els based on deeper data, object subtracted, edited to remove 
artefacts, background-subtracted) on the GUViCS web sit^jand 
make them publicly available. We also make available the com¬ 
plementary maps on the same 1 arcmin pixel scale prepared and 
used for our analysis. 
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